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■^The  equilibrium  behavior  of  Jackson  queueing  networks  (Poisson  arrivals, 
exponential  servers  and  Bernoulli  switches)  has  recently  been  investigated  in 
some  detail.  In  particular,  it  was  found  that  in  equilibrium,  the  traffic 
processes  on  the  so-called  exit  arcs  of  a Jackson  network  with  single  server 
nodes  constitute  Poisson  processes.  This  result  m.-Tv  be  viewed  as  an  exten- 
sion of  Burke's  Theorem  from  single  queues  to  networks  ol  queuos. 

A conjecture  ma d e by  Burke  and  others  contends  that  the  traffic  processes 
on  nonexit  arcs  cannot  be  Poisson  in  equilibrium.  This  paper  proves  this  


FORM 
I JAN  73 


EDITION  OF  < NOV  84  IS  OBSOLETE 


Unclassified 


Une lass i f ied 


w W " J Jl 


i hi  j P AOE'^/'pn  Data  Entered) 


conjecture  to  be  true  for  a variety  of  Jackson  networks  with  single  server 
nodes.  Subsequently,  a number  of  characterizations  of  the  equilibrium 
traffic  streams  on  the  arcs  of  open  Jackson  networks  emerge,  whereby 
stochastic  properties  of  traffic  streams  are  shown  to  be  equivalent  to  a 
simple  graph-theoretic  property  of  the  underlying  arcs.  These  results 
then  help  to  identify  some  inherent  limitations  on  the  feasibility  of 
equilibrium  decompositions  of  Jackson  networks. 
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Abstrac  t 


The  equilibrium  behavior  of  Jackson  queueing  networks  (Poisson 
arrivals,  exponential  servers  and  Bernoulli  switches)  has  recentlv  been 
investigated  in  some  detail.  In  particular,  it  was  found  that  in  equil- 
ibrium, the  traffic  processes  on  the  so-called  exit  arcs  of  a Jackson 
network  with  single  server  nodes  constitute  Poisson  processes.  This 
result  may  be  viewed  as  an  extension  of  Burke's  Theorem  from  single 
queues  to  networks  of  queues. 

A conjecture  made  by  Burke  and  others  contends  that  the  traffic 
processes  on  nonexit  arcs  cannot  be  Poisson  in  equilibrium.  This  paper 
proves  this  conjecture  to  be  true  for  a variety  of  Jackson  networks 
with  single  server  nodes.  Subsequently,  a number  of  characterizations 
of  the  equilibrium  traffic  streams  on  the  arcs  of  open  Jackson  networks 
emerge,  whereby  stochastic  properties  of  traffic  streams  are  shown  to 
be  equivalent  to  a simple  graph- theoretic  property  of  the  underlying 
arcs.  These  results  then  help  to  identify  some  inherent  limitations  on 
the  feasibility  of  equilibrium  decompositions  of  Jackson  networks. 

/ 


1.  Introduction. 


P.J.  Burke  [4]  and  E.  Reich  [14]  showed  independently  that  the 
departure  process  from  a M/M/s  queue  in  equilibrium  is  a time  homogen- 
eous Poisson  process  (with  same  rate  as  the  arrival  process) . Analo- 
gous results  have  been  recently  obtained  by  F.P.  Kelly  for  a large  class 
of  queueing  networks  with  exponential  servers  ([11],  [12]),  and  even  for 
certain  ones  with  general  independent  servers  [1].  These  results  made 
elegant  use  of  reversibility,  first  employed  by  Reich  in  [14].  Similar 
results  were  attained  more  directly  ([3];  [13]  Sec.  4.7)  for  the  class 
of  Jackson  queueing  network.  These  networks  are  termed  after 
R.R.P.  Jackson  [10]  and  J.R.  Jackson  [9];  they  are  described  below. 

A Jackson  network  consists  of  m service  stations.  Each  station  i 
houses  a finite  number  of  parallel  identical  independent  servers  that 
provide  exponential  service  times  with  rate  .3  ^ . Customers  can  arrive 
at  a station  i either  endogenously  (from  other  nodes)  or  exogenously 
(from  outside  the  network);  in  the  latter  case  they  arrive  from  eacli 
source  according  to  a Poisson  process  with  intensity  When  the 

servers  are  busy,  customers  are  delayed  in  a FCFS  (first  come  first 
serve)  waiting  line  of  infinite  capacity.  On  service  completion  each 
customer  enters  a Bernoulli  decomposition  switch  where  his  next  desti- 
nation is  determined  according  to  a multinomial  Bernoulli  trial.  Farh 
customer  may  then  he  routed  to  station  j with  probability  p..,  or  leave 

t 

tin'  system  altogether  with  probability  q I - ) p...  All  exogenous 

1 M 11 

arrival  processes,  service  times  and  switching  decisions  are  mutually 
i n dependent . 

In  this  paper,  we  shall  he  solely  concerned  with  Jackson  networks 
with  single  server  stations.  Surh  networks  are  compactly  specified  bv 
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1 

J 


a quadruple  JN  = (M,a,rJ,P)  where  M = {l,2,...,mJ  is  the  station  set, 

a = (a, , . . . ,u  ) is  the  vector  of  arrival  rates,  o = (a,  , . . . ,<>  ) is  the 
1 m 1 m 

vector  of  service  rates,  and  P = [ p . . ] is  a b x in  substochastic  matrix 

ij 

which  specifies  the  switching  probabilities  among  stations.  The  dummy 

A 

station  0 will  denote  the  network  sink,  with  the  convention  p^  = 1. 

Since  a network  configuration  is  routinely  pictured  as  a directed 

graph,  we  shall  find  it  useful  to  couch  much  of  the  impending  discussion 

in  graph-theoret ic  terms  (c.f.  [151).  Accordingly,  service  stations 

will  be  termed  noth;::,  while  the  term  arc  will  refer  to  node  pairs  (i,  j) 

such  that  p..  > 0.  A path  is  any  sequence  of  nodes  (k  ,k  , ...,k  ) such 
ij  1 2 n 

that  p , >0  for  every  1 < l<  n.  The  lenqth  of  a path  (k  ,k  ,...,k  ) 

H.+l  ~ ' 1 2 n 

is  n-1,  i.e.  the  number  of  adjacent  arcs  specified  by  the  path.  If 

u.  > 0 then  the  network  is  said  to  have  an  inlet  to  node  i.  Likewise, 
i 

an  arc  (i,0)  is  said  to  be  an  outlet  of  the  network  from  node  i,  if 
> 0. 

The  accessibility  relation  is  defined  on  MU{0)  as  follows:  j is 

accessible  from  i (notation  i j ) if  i = i or  there  is  a path  (i,k  ,...,k  , 

I n 

We  say  that  i and  j communicate  if  i j and  j - v i.  The  communication 
relation  is  an  equivalence  relation  on  MUfOl;  the  equivalence  classes 
of  the  resultant  partition  will  be  termed  here  cnmpenrnls  of  the  network. 

An  elementary  classification  of  nodes  and  Jackson  networks  now 
follows.  A node  i is  <>prn  if  i~-»0;  it  is  il  \sf~ttt.  A Jackson 

network  is  i<y;  or  •/  •.•<  / according  as  all  its  nodes  .ire  open  or  closed; 
otherwise,  it  is  mi.r>:,l.  A Jackson  network  is  ,■/</<  e m if  it  has  no 
inlets  (see  Sec.  4.5  in  fill  for  more  details).  A fundamental  classi- 
fication of  arcs  is  given  by 
he  I i n i t i on  I . 1 (cl.  I 1 | ; | I 1 | Sec . 4.7) 

An  arc  (i,j)  is  an  exit  arc  if  j.)L»i.  Otherwise  (i,j)  is  a nonexit 
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Thus,  an  exit  arc  is  located  between  distinct  components,  while  nonexit 

arcs  reside  within  components.  It  also  follows  that  a nonexit  arc  (i,j) 

has  a cycle,  viz.  there  is  a path  of  the  form  (i,j,k,...,k  ,,i).  Con- 

1 n-1 

sequently,  (i,j)  has  a cycle  of  minimal  length  L=L(i,j);  such  cycles 
will  be  termed  minimal  cycles  (of(i,j)).  In  particular,  when  L=l,  (i,i) 
will  be  called  a feedback  arc. 

Let  {K  (t)  be  the  traffic  (counting)  process  on  arc  (i,  j), 

viz.  K^.(t)  it  the  total  number  of  customers  that  passed  on  arc  (i,j) 
in  the  time  interval  (0,t].  The  salient  feature  of  traffic  processes 
on  exit  arcs  (i,j)  is  that  in  equilibrium  {k„(t)lt>Q  is  a Poisson  pro- 
cess. (See  [3];  [13]  Sec.  4.7).  A stronger  result  will  be  stated  later 
on  in  'theorem  2.1. 

The  study  of  traffic  processes  is  essential  to  the  study  of  network 

decompositions.  it  was  conjectured  in  [13]  (Conjecture  4.7.2)  that  if 

(i,j)  is  a nonexit  arc  (excluding  feedback  arcs  will)  p..=]),  then 

[ K . j ( t ) f t ()  in  not  a Poisson  process.  A similar  conjecture  was  previousL 

asserted  by  P..J.  Ilurke  [5],  where  a heuristic  argument  to  that  effect  is 

also  outlined.  Preliminary  results  supporting  the  non-Poisson  conjecture 

were  obtained  in  [13]  (see  Sec.  4.7). 

In  this  paper,  we  shall  prove  a refined  version  of  the  atoresaid 

conjecture  as  regards  feedback  arcs  in  arbitrary  Jackson  networks;  we 
shall  prove  it  to  be  true  for  any  open  Jackson  network  in  equilibrium. 

In  the  latter  case,  this  will  yield  a Poisson  characterization  of  the 
traffic  processes  based  on  the  exit  property  of  the  underlying  arcs. 

A number  ot  other  characterizations  will  be  seen  to  emerge  from  this 
result;  these  results  suggest  general  principles  that  impose  inherent 
limitations  on  network  decompositions.  In  addition,  we  shall  seek  to 
shed  some  light  on  the  nature  of  the  deviation  from  the  Poisson,  and  to 
point  out  the  interplay  between  the  graph-theoretical  aspects  and  the 


statistical  aspects  of  Jackson  networks. 
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Review  of  So  ini'  Known  Results 

In  this  section  we  review  some  known  results  that  pertain  to  Jackson 


networks; the  term  Jackson  networks  will  henceforth  refer  to  networks 
with  single  server  nodes. 


Let  JN=(M,a,a,P)  specify  such  a network.  The  associated  bntiffi 


i/uaLi<>n  is 


(2.1) 


t = u + t P 


in  the  unknowns  6 =(d  , . . . ,6^) . A solution  <S^0  for  (2.1)  is  called  a 

traffic  solution.  An  open  Jackson  network  always  has  a unique  traffic 

solution;  an  autonomous  closed  one  has  infinitely  many  solutions  (see 

[2], and  [13]  Sec.  4.4  for  more  details). 

The  state  process  of  the  network  is  denoted  (Q(t))t>Q  where 

Q(t)=(Q, (t) , . . . ,Q  (t))  is  the  vector  of  customer  totals  at  each  node  at 
1 m 

time  t.  (Q(t) } q is  a Markov  process  ([2];  [13],  Theorem  4.2.2)  whose 

state  space  ranges  over  all  m-tuples  v=(n,,...,n  ) of  non-negative 

i m 

integers. 

For  each  V=(n,,...,n  ) denote  for  brevity  P [Q. (t)=n, , . . . ,Q  (l)=n  ] P ( ), 
1 m 11  m m t 

and  let  £ . denote  the  m-dimensional  unit  vector  with  1 in  the  i-th  coor- 
l 

dinate.  The  birth-and-death  equations  of  (Q(t)}  q are  (cf.  [9]) 

(2.2)  ~~  P (v)  = 

f)  t t 

m 

) P , ( V>— i .)<«.. 

1=1  L 1 1 


+-  ) P.  (v+f  . >J  .q  . 

j=l  1 1)1 


4 > > P (V-i  .+<  . )0  .p  . . 

f ' I IP 


m m 


- P,  (v)  ( ).  '1.  F ) o .q  . ( (n  . )+  y } o . p . . | (n  . ) ) 

' 1=1  ' j=l  11  1 |=1  


V=(n,  , . . . ,n  ) (),  t ’0  , 
m - - 


I ip--;' 
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where 


f ( n ) = 


0,  if  n=0 


L,  if  n>0 


A Jackson  network  is  said  to  be  in  equilibrium  if  JQ(t)  ^ is  in 

3 ~~ 

steady  state;  equivalently,^—  P (v)~0  V V=(n  , ...,n  )>0. 

d t t 1 m - 

A Jackson  network  in  equilibrium  always  has  a traffic  solution 

A m 

6 = (E [ D ( t , t+1 ]],..., E [D  (t,t+l] ]), where  D.(t)  = V K..(t)  and 
1 m l i 


on 


A J=0  1J 

D^(t,u]  = D^(u)-D^(t)  ([2];  [13],  Theorem  4.5.4);  hence, 6 may  depend 

the  initial  condition  P (v) , v>0.  The  relation  between  the  Traffic 

o - 

Equation  and  network  equilibrium  is  investigated  in  detail  in  Sec.  4.5 
of  [13];  in  particular,  if  the  network  is  not  open,  it  is  necessary  in 
equilibrium  for  each  closed  part  to  be  autonomous  (ibid.). 

For  mixed  networks  in  equilibrium,  the  latter  fact  enables  us  to 
deal  separately  with  each  open  part  and  each  closed  autonomous  communi- 
cating part. 

For  open  Jackson  networks,  the  equilibrium  state  distribution  was 
shown  by  J.R.  Jackson  [9]  to  be 

m n . 

(2.3)  P[Q.  ( t)=n  Q (t)=n  ] = II(1-P  IP.1 

1 1 m m . , l l 

i = l 


A 'J 

under  the  sufticient  condition  1,  i 1 M . The  aforesaid  condition 

i . 

is  also  necessary  for  equilibrium  ([2];  [13]  Theorem  4.5.5).  Similarly, 
for  autonomous  closed  networks  in  equilibrium  such  that  the  number  of 
total  customers  in  the  network  equals  n with  probability  1,  W.J.  Gordon 


.and  G.F.  Newell  |B|  showed  that  in  equilibrium 


(2.4) 


P[Q.  (t)=n  ()  ( t )=n 

11  mm 


SM  "I  g(r») 


m n 

!l  t J 

I 


i = l 


where  t , = * for  some  traffic  solution  , and  g(n)  is  a normalization 

l <5  . 

1 

constant.  Notice  that  the  state  space  is  restricted  to  >-(n,,...,n  ) 

m 

m 

such  that  ) n.=n.  More  generally,  if  an  arbitrary  distribution  ot  S 
i = 1 1 
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is  known,  then  in  equilibrium 

(2.5)  I'[Q.  (t)=ni  , . . . ,Q  (t)=n  | = 

J 1 mm 

m m 

P[Q,  (t)=n  , . . . ,Q  (t)=n  | S = ) / n.]. 

1 1 m m 1 M . i M i 

1=1  i=l 

Let  (a, b)  be  an  arbitrary  arc.  Consider  the  process 

{(Q(t);  Kgk  ^ t>o  l:^e  state  augmented  by  the  traffic  count  on  arc 

(a,b).  { (Q ( t ) ; K^Ct)))  >q  is  a Markov  process  with  a countable  state 

space  of  the  form  (v;k)>0  ([3];  [13]  Theorem  4.7.1).  In  writing  the 

corresponding  bir th-and-death  equations, we  must  distinguish  between 

two  cases.  Here, the  C . are  the  (mfl)-dimensional  unit  vectors  and 

P (v;k)  = P[Q  ( t)=n  , . . . ,Q  (t)=n  ,K  ( t ) =k  ] , for  any  ( .;k)  = (n  > ,k)  0. 

t I 1 m m ab  1 m - 

Case  1 (feedback  arcs) : a = b 

m 

(2.6)  y-  P (v;k)  = l P ((v;k)-e.)cx. 

i=l  11 

m 


+ I P ((v;k)+e.)a .q. 

j=l  t 31  3 

m m 

+ 11  pt ( (v; k)-&i+e . )a . p 

1=1  1=1  J 3 

j/i 


+ Pt((^;k)-c^l)aap 


m m 


a a 

m m 


- P ( V ; k ) [ l a + |o  q f(n  )+  1 1 a p f (n  )4  j p f(n  )], 

i=l  j=l  1 1 1 i=]  j=l  J ]1  1 J 33  3 

jti 

(v»;k)  >0,  t >0. 


Cast*  2 (nonfeedbuck  arcs):  n/b 


^ • asas^r  3S3BES 


+ Pr((v;k)-e  +e  -e  )cr  p 

t b a m+1  a ab 


m m 


> (v;k)[  l a + l a q P(n  )+  l l a p B(n.)], 
i=l  j=l  3 3 3 i=l  j=l  J J1  J 

j*i 


(v;k) >0,  t>0 . 


The  initial  conditions  in  each  case  are 

(2.8)  P (v;k)=  ! P (v) , if  k=0 
()  jo 


\ 0 , otherwise. 

Finally,  we  cite  two  results  pertaining  Lo  equilibrium  traffic 


processes . 


Theorem  2.1 


Let  JN-(M,a,a,P)  be  a Jackson  network  in  equilibrium.  Let 
C ^-1  * j i ) > • • • > ( i-n » J n ) be  a set  of  arcs  such  that  j rfo  * for  all  KSl<n, 


1 • kT  n . 


Then  the  traffic  processes  {K.  . (t)}  , 


. • (t)  \ , 1-  r n are  mutually 

l t^o  - - ’ 

r r 


independent  Poisson  processes. 


See  [j]  or  [131,  Theorem  4.7.3. 


Theorem  2 . 2 

In  equilibrium,  each  component  C of  the  network  constitutes  itself 
a Jackson  network. 


See  [2  | or  [13],  Theorem  4.7.4. 
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3.  Some  Implications  of  a Poisson  Traffic  Process 

A m ni  • 

Let  <f>  (z  , ...,z  ) = V P (n  , ...,n  )(  H z ) be  the 

tl  m f u v.  tl  m . , 1 

(n  , ...,n  )>0  i-l 

1 m - 

generating  function  of  the  distribution  of  Q(t).  Likewise,  the  gener- 
ating function  of  the  distribution  of  (Q(t)  jK^Ct))  is  denoted 

A 3 m ni  k 

'!>  (z, z.y)  = y P (ni  , . . . ,n  ,k)  ( II  z ) y . 

1 1 m (n,  , . . . , n , k)  >0  t 1 m 1=1  1 

The  generating  functions  above  and  in  the  sequel  are  defined  for  I 0, 


\y.  I*  l (i»M),;md  |y  * I, unless  olh(*rwise  specified.  We  shall,  « » I course, 
i - - 

use  the  convention  0°  = 1.  Notice  that  <p  (z , 1 z ,y) 

i 1 l m 

=4>t(zi> • • • *zi_i» zi+i> • • • >zm»y)  where  1 designates  the  substitution  of 
1 for  in  the  argument  list.  This  is  so,  because  setting  z.=l 
results  in  the  generating  function  of  the  respective  marginal  distri- 
bution by  Abel's  Theorem  (cf.  [9],  1.2).  In  general,  substituting  a con- 
stant c for  z.  in  the  argument  list  will  be  denoted  ip  (z, ,... ,c .,..., z ,y). 

i t 1 i m 

The  probability  trajectories  P (n,...,n  ,k)  have  derivatives  of  every 

t 1 m J 

order  in  t.  Moreover,  every  countable  sum  thereof  is  uniformly  con- 
vergent on  each  compact  interval  of  [0,no).  This  fact  will  justify  all 
termwise  operations  to  be  performed  in  the  sequel,  such  as  termwise 

passage  to  limits,  integration ,dif ferentiation,  etc.  (cf.  [9],  1.1,  1 . 7 » . 


We  now  proceed  to  write  the  generating  function  version  of  (2.6) 
and  (2.7).  By  adding  and  subtracting  the  term  P ( (v  ;k)-r  ^ )o  ^p  ^ in 

Case  2,  we  may  combine  the  two  cases  into  one  case  as  follows 


(i.D  vv-’vy) = y vv---’vy)(ii(vi)  + 


i=i 


+ } , (lf't (V  • • • ’ 1 0 j ’ • • • * Vy) ),T  jq  j (z7  _1  > 


m m 


1 >.  ). 

i i j=i 

i/i 


,z  , y ) — '|>  (z  ,...,0 z ,y))i.|i..(  -I) 

in  1 I I m | | i z . 


+ (^t(zl’'"*Vy)-\W*'--’0a’--*’Vy))aaPab(—  “1} 

a 

Zb 

- (<>t(z1, . . • » zm» y)_^t ( Z1  ’ • • • ,0a’  ■ * ■ ,Zm,y))aaPab(z  1)  ' 


Note  that  the  singularities  due  to  z.=0  in  the  respective  denomi- 
nators are  all  removable  via  continuous  extensions. 

Next  set  z^=l,  irM,  on  both  sides  of  Equation  (3.1).  We  obtain 

( 3 - 2 ) ^ lPt(y)  = (<Pt(y)-<Pt(Oa,y))oapab(y-l) . 

In  the  sequel,  {B^(t)}t>Q  will  denote  the  state  indicator  process  of 


node  i,  where 


b.(t)  = 


0,  if  Q . ( t)=0 


^1,  if  Q.(t)>0. 


Lemma  3 . 1 


(3.3)  E[Kab(t)]=OaPab  • / Pf Ba(u)=J ]du,  t>0. 


Proof 

Equating  coefficients  in  Eq.  (3.2)  resul ts , af ter  some  manipulation, 
in  the  system  of  equations 

fc  P[Kab<t)-kl* 

(P[Ba(t)=l,Kab(t)=k-l]-  P[Ba(t)=l,Kab(t)=k])  Japab  , 
k=0, 1 

For  each  fixed  n 1 integrate  and  nn  holli  sides  above  over  k n. 


I II  view  of  (2.8)  we  ;>c  I 

l 

,>[Kab(t)^n]=  / P(  I$a(u)=l,K  ib(u)=n-l  | o^p  du  , n I. 
o 

Since  K , (t)  is  a non-negative  integer  valued  random  variable 
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E[Kh(t)l-X  P[K,(t)>n]= 


).  / |MBa(»)-1,Kab(,i)=i.-l]«i|P{|b«|„= 

n=l  o 

t 

Vab  / PfBa(u)=l]du 

o 

by  interchanging  summation  and  integration.  n 

Corollary  3.1 

E[Ka^(t) ]=At,  t>0,  for  some  A>0  iff  {Ba(t)}t>g  is  in  steady  state. 
In  particular,  if  (Kala ( t ) } ^ >q  is  a time  homogenous  Poisson  process, 
then  {B  ( t)}  „ is  necessarily  in  steady  state. 

a t [ | 

We  can  now  characterize  Poisson-distributed  traffic  processes  in 
Jackson  networks. 

Theorem  '3.  1 

K , (t)  is  Poisson  distributed  for  every  t>0  iff  13  (t)  and  K , (t) 
ab  - a ab 

are  independent  for  each  fixed  t>0. 


Proof 


(=*•)  Suppose  that  K^^(t)  is  Poisson  distributed.  Thus 


.and  consequently  from  (3.3) 


;i(  'l>t(y)=,l't(y),’aP.lb  P[H;i(t)=l l(y-l)  . 

Substituting  the  above  on  the  left  side  ol  (3.2)  ami  s i mp  I i I y i ng 


y i e I (Is 


t(y)  P[Ba(t)=l  ]=<|-t(y)-l' t(«;1, y),  |y|-  1- 


Equating  coefficients  on  both  sides  above  results  in 

P[K  , ( t )=k  | * I* | B ( t ) =1  ]=  P|K  . (t)=k  |-  P [Q  ( t )=0, K , ( t )=k  | 


p[Knb(t)=k,Ba(t)=l |,  k=0 , 1 , 
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The  requisite  independence  now  follows, since  Bg(t)  is  a zero-one  random 
variable. 

(*)  Suppose  Ba(t)  and  Kab(t)  are  independent  for  every  fixed  t. 

Then  <J>  (0  , y)=<}>  (0  ) <J)  (y)  ,and  Equation  (3.2)  becomes 
L 3 t a C 

4>t(y)=(^t(y)-4>t(oa)<0t(y))aapab(y-1). 

This  can  equivalently  be  written  as 

ft  Vy  } = Vy } ‘ P 1 Ba  ( c } =1 1 °apab  (y_1 } * 

and  the  unique  solution  is 


4>t(y)« 


eP[Ba(t)=l]aaPab(y-l)> 


Finally,  the  above  is  recognized  for  each  fixed  t as  the  generating 
function  of  a Poisson  distribution  with  respective  parameters 


P[Ba(t)=1]aaPab- 
Corollary  3.2 

If  {Kab(t))t>Q  is  a time  homogenous  Poisson  process, then  in 
part  i cilia r 


""n  P[Q  (t)=0,Kab(t)=0J=  P[Qa(t)=0]e"MKab(,)  |l(-E[Kab 

ft  t 

The  following  theorem  characterizes  time  homogenous  Poisson 
in  Jackson  networks. 

Theorem  3.2 


Dl) 

1 1 

traffic 


(K  ,(t)i  is  a time  homogenous  Poisson  process  iff 
an  t>0 

a)  lK,1j)(t)}j  j)  Is  a renewal  process 
and 


fli./Dl 


t 0 


l>) 


Is  In  steady  slate. 
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I 

Proof 

(=5-)  If  {K  , ( t ) } _ n is  a Poisson  process, then  it  is  well-known  that 
ab  t'MJ 

a)  holds,  and  b)  follows  from  Corollary  3.1. 

(*=)  Suppose  a)  and  b)  hold.  From  Corollary  3.1  it  follows  that 

R(t)  4 E[Kab<t)l-OaPab  P[Ba(0)-l]t 

where  R(t)  is  the  renewal  function  of  (K  , (t)}  But  the  only 

ab  t>U 

renewal  process  with  R(t)=At  for  some  A>0  is  the  Poisson  process 
([6]  p.308).  □ 

It  should  be  pointed  out  that  some  aspects  of  the  results  obtained 
so  far  are  strongly  reminiscent  of  analogous  ones  attained  in  [7]  in 
the  context  of  a single  M/G/l  queue. 

The  reader  is  reminded,  however,  that  our  results  apply  to  any 
Jackson  network  with  arbitrary  topology  (open, closed  or  mixed).  With 
the  aid  of  these  results  we  shall  proceed  to  treat,  in  the  sequel, 
traffic  processes  on  nonexit  arcs.  The  discussion  will  be  restricted 
to  arcs  (a,b)  which  are  traffic-nontrivial  in  the  sense  that 
E[Kab(t)]?0,  t>0,  or  equivalently  Pab' p [ Ba ( t )=1 ]70,  t>0. 
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4 . Traffic  on  Feedback  Arcs 

Throughout  this  section,  (a, a)  will  designate  a traffic-nontrivial 
feedback  arc. 

Before  proving  the  main  theorem  we  need. 

Lemma  4 ■ 1 

If  P[B  (t)=in  t>0,  then  p =]  . 

3 — 33 

Proo  f 

Setting  z =1  for  all  i'M-faf  in  the  generating  function  version  of 
(2.2)  yields  (cf.  12]), 


§7  VZa)  = 


(VZa)aa+.^  t(za)^t(za’°j))ajPja 

+ (*t(za)-»t(Oa))oa(l-Paa)(f-l). 


)<v« 


The  inverse  transformation  of  the  above  to  the  time  domain  gives 
after  some  manipulation 

(4.1)  P[Qa(t)=n]  = 

(P [ Q_ ( t ) =n-l ]-  P[Q  (t)=n])a 

3 3 3 

111 

+ y ( P[Q  (t)=n-l,Q.(t)>01-  P[Q  (t)=n,Q.(t)>01)  7 p. 

a J a J J Ja 


+ 0,[(J.i(t)=n+l  I-  P[Qa(t)=nl+  P[Q;l(t)=0|.^(i>v  (1-p  ( ,), 

n=0, 1,2,... 


wliere  il  is  Kronecker's  delta, 
no 

Now, suppose  I’fB  (t)  = l I 1,  t O.but  p • I.  We  show  hy  indiiel  ion  on 
n=0 , 1 , . . . t ha t 


(4.2)  P [Q  (t)=k]-0,  t>0,for  all  k< n 

3 “ ~ 

which  will  contradict  P[B  (t)=l]  1,  t>0.  For  n=0,  (4.2)  immediately 
follows  from  the  supposition. 

Assume  that  (4.2)  holds  for  n>0.  In  view  of  the  induction 
hypothesis  we  obtain  from  (4.1) 

0=  P[Qa(t)=n+l}Ja(l-paa) 

which  establishes  the  induction  step  since  0 (1-p  )?^0  by  the  supposi- 

3 33 


Theorem  4.1 

Let  JN= (M,a,a ,P)  specify  an  arbitrary  Jackson  network.  Then 
{K  (t)>  n is  a time  homogenous  Poisson  process  iff 

33  t 

a)  P[B  (0)=1]=1 

a 


b)  p =1 
1 na 


Proof 


For  a=b  send  t >0+  in  Equation  (3.1)  thus  obtaining 


(4.3) 


• it  V*! Vy)  = 


I <Do(z1,...,Zni,y)ai(Zi-1) 

1 = 1 


m , 

+ I (Vv-v^'V2! °j vy))Vj(z: _1) 

j=l  j 

mm  7, , 

+ l \ (<t>  (z, z ,y)-<f>  (z o ...,Z  ,y)0.p..(  1 -1) 

i = l : = l ° 1 m o 1 j m I M z. 

‘ 

f (<()  (/  ,...,/.  ,y)-'|i  (z  ,y))i  p (y-l). 

o I in  o I a m a aa 

In  view  of  (2.8), the  transformed  version  of  the  initial  eonditii 


Notice  that  in  view  of  l.etnma  4.1,  a)  and  b)  above  are  equivalent  to 
P f B ( ( t )=.l  1 1,  t_0. 
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(4.4)  <p  (z  ,...,z  ,y)-<J>  (z  , . . . ,z  ) . 

o 1 m o 1 m 

A comparison  of  (2.7)  in  Case  1 with  (2.2)  reveals  that  in  view  of 

(4.4) ,  Equation  (4.3)  reduces  to 


K-  4>  (z  > • • • .Z  ,y)  = 
d t o 1 m 


K—<p(z,...,z  ) + (<p  (z  ,...,z  )-<p  (z  , ■ ..,  0 , . . . , z ))o  p (y-1 ) 
otol  m ol  mol  a ma'aa7 


Setting  in  (4.5)  z^l,  ifM-fa},  yields 


(4.6)  - (J>  (z  ,y)  = r—  $ (z  ) + (<k  (z  )-<]>  (0  ))n  p (y-1) 

> t o a a t o a o a o a a aa  J 


(=*•)  Suppose  IK  (l)l  is  a time,  homogenous  Poisson  process.  On 
aa 

setting  z=0,  y=0  in  (4.6)  we  observe  that 
a 

(4.7)  ~ P [ B (0)=0, K (0)=0 ]=  P[Q  (0)=0]+0=0 

* u d aa  a t a 

since  P[Q  (t)=0]=-0,  t>0,by  Corollary  3.1.  But  from  Corollary  3.2 
* t a — 


(4.8)  ' P[B  (0)=0,K  (0)=0]=  P [Q  (0) =0 ] (-E [K  (1)  ] ) . 

» i a aa  d aa 

From  (4.7)  and  (4.8)  we  deduce  p[Q  (0)=0]=0,as  E[K  ( 1 ) ] >0  by  traffic- 

a aa 

nontriviality,  so  that  Condition  a)  is  established.  Furthermore,  by 
Corn! 1 a ry  3.1 


P[Ba(t)  = l|  1,  Id). 

Condition  b)  now  follows  from  Lemma  4.1. 

(*)  If  Conditions  a)  and  b)  hold,  then  lKa-j(t)}  ()  is  a time 
homogenous  Poisson  process ,s ince  the  interdeparture  intervals  are 
mutually  independent  and  exponent  ially  distributed  random  variables 
with  same  parameter  o . 

I Observe  that  (4.5)  is  valid  even  when  we  deal  with  a singleton  auton- 
omous closed  network  (one  with  M=  la  I and  pna  = l ) . In  this  case,  all 
terms  on  the  right  side  ol  (4.  !)  except  for  the  last  one  are 
empty  and  therefore  evaluate  to  zero.  However,  it  can  be  directly 

verified  that  ^ i|(  (z  ) 0,  t>0,  as  required. 


i 


I 
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5.  Traffic  on  Non  feedback  Nonexit  Arcs 

Let  (a,b)  be  a traffic-nontrivial  nonfeedback  nonexit  arc;  thus 

a^b ,and  the  minimal  cycle  length  (see  Sec.  1)  satisfies  L(a,b)=L>l.  Our 

working  conjecture  is  that  ' ( t) } Q can  never  be  a Poisson  process. 

The  proof  strategy  would  be  analogous  to  the  one  in  the  previous  section, 

rramely  to  demonstrate  an  inconsistency  with  Theorem  3.1.  This  in  turn 

would  again  enable  us  to  deduce  that  {K  , (t)}  . „ is  neither  Poisson 

ab  t ->0 

distributed  nor  is  it  a renewal  process.  A more  extensive  conjecture, 
consistent  with  the  above,  is  put  forth  in  the  following 
Conjecture  5.1 


p[Qa(0)=0,Kab(0)=0|= 


P[Qa(0)=0](-<Japab  P[Ba(0)  = l])L 

+ l d(a,j  ) 

C . (a.b)  1 L 

mm 

where  the  sum  is  taken  over  all  (a,  j j T ) in  the  set  C . (a,b)  of  ill 

1 L nun 

minimal  cycles  of  (a.b),  and  the  d (a , j , . . . , j ) are  nonvanishing  terms 
all  having  the  same  sign. 


Prooving  this  conjecture  would  necessitate  the  computation  ot 

I. 

, if  ( 7 ; ,v).  We  now  outline  how  this  may  be  done. 

(.  o 1 m 

a t 

Let  f>  ^ be  the  set  of  generating  functions  (z- transforms) 

y(z,  .....  z ,y)  of  sequences  {J'(l0  I ,,%ri  .where  p=(n  ,,..., n ,k)  ranges  over 
I m |i  u I m 

all  (nri-l)-tuples  of  non-negative  integers,  and  such  that  /(it)1'". 

1*  jo 

Next,  associate  witli  .IN-  (M,  't,n  , P)  ami  arc  (a,b)  an  operator 

„ . . ( m+ 1 ) , ,(m+ 1 ) . . 

I):  (.  mi  def  ined  hv 


(5.  I ) l)|  y(z, /.  , v)  I 

m 


Z , V ) o ( z - I ) 

111  I | 


1 
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m 1 
+ l (y(z  ,y)-y(z  , . . . ,0. , . . . ,z  ,y))o  .q . (—  -1) 

1 m 1 i m l i z . 

j=l  j i j j 

mm  7 . 

+ I I (y(z  ,...,z  ,y)-Y(z  ... ,0 z ,y)«.p..(-‘  -1) 

i=1  j=l  1 m 1 1 tn  , z. 

yti 

+ (Y(z1,...,Zm,y)-Y(z1,...,0a,...,zm,y))Oapab  (y-1)  , 

a 

| z . |<1  (ieM)  , | y I < 1 . 

l - - 

We  make  the  following  observations 


1. 

n ‘M2!  » ‘ • 

at  1 

• , z , y ) =D°  [ <f> 
m 

(z  ,...,z  ,y)]=Dt‘ 
t ! m 9 tn 

2. 

lim  D[<}>  (z..  , 
t*0+  c 1 

...,z  ,y)  ]=D[lim  <(|  (z  , . ..  ,z  ,y)] 
m t->o+  1 m 

whence 

3. 

Dn[<()  (z  , . . . 
o 1 

,y)  ]=  — - 

at 

Vzi vy) : 

and  finally 

4.  D is  linear  viz. 

D[clVC2Y2,=ClD[Yl  ]+C2D^Y2  ^ 

for  any  ‘Yi»Y2e®  and  scalars  c^,C2* 

Observation  1. reflects  the  fact  that  on  differentiating  both  sides 

3n 

of  (3.1)  n-1  times,  <b  (z,  , . . . , z ,y)  can  be  obtained  recursively 

i -n  „n  t 1 m 

^ n-1  o t 

from  — — - <p (z, , . . . , z ,y)  by  setting  the  latter  on  the  right  side 
n-1  t 1 m - 

i t ^n-i 

of  (3.1)  (or  equivalently  by  applying  D to  r-  <P  (z,  , . . . , z ,y)). 

-i  n-1  t 1 m 

o t 

In  part icular, this  is  true  for  t=0  due  to  the  continuity  of  I)  in  the 

sense  of  Observation  2.  Thus,  the  problem  of  computing 
( I. 

, 1 .y)  reduces  to  Lite  combinatorial  problem  of  calculating 

,il  ' ° 1 m 

I)1' fill  ( z z , v | . A major  reduction  in  the  combinatorics  involved 

o l m 

can  be  attained  bv  collecting  terms— much  like  in  the  previous  section  — 
and  completing  them,  if  necessary,  to  expressions  which  vanish  under 
the  appropriate  conditions.  Thus,  sending  t >()+  in  (3.1)  lor  a^b  yields 
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(5.2) 


I*  (z  > . . . ,z  , y ) | 
o 1 m 


^ b 

FT  <J>  (z, ) + (A  (z  , . . . ,z  ,y)-<b(z  , . . . ,0  . . ,z  ,y))  ) p ’ (y-1) 

d t o 1 m ol  m ol  a m a an  z y 

a 

g 

due  to  (2.8).  Notice  that  r — 4>  (z,,...,z  )=0  in  equilibrium. 

o t o 1 m 

Lemma  5 ■ 1 

For  {Kab(t) } t>()  to  be  a Poisson  process,  it  is  necessary  that  for 

every  n'0,  --  P[Q  (0)=n,K  , (0)=0]  = ~ P [Q  (0)  =n  ] -P  [ Q (0)=n+l  | T p . 

*“  o l a an  i)  i d a a at) 

P roof 

Set  z.=l  f°r  all  ij^a  in  (5. 2),  thus  obtaining 

'v  I'  (*  >y)  = A 'I'. (z  )+(<l*  (z  (0  )>!  P , 1 (y-I)  . 

It  o a It  o a o a o a a ab  z 

a 

The  Lemma  follows  by  setting  above  z_^=0  and  y=0. 
in  view  of  Corollaries  "3.1  and  3.2  we  have 


If  is  a time  homogenous  Poisson  process, then  necessarily 


P[Q  (0)  = 1 ]=  P [Q  (0)=0 ] * P [Q  (0) >0 ] 

a a a 


P [Q  (0)>1]=  P [Q  (0) >0 ] *P  [Q  (0)>0] 
a a a 


.'orollary  5.1  considerably  restricts  the  feasibility  of  Poisson 
traffic  processes.  An  inspection  of  Equations  (2.4)  and  (2.3)  discloses 
that  even  in  equilibrium  a closed  network  floes  not,  in  general,  sat  islv 

the  necessary  conditions  of  Corollary  5.1.  However,  I rum  (2.  1),  m open 

Jackson  network  in  equilibrium  always  does.  In  view  ol  the  feasibility 
of  decomposing  the  latter  networks  according  to  Theorem  2.2,  we  shall 
now  proceed  to  show  that  Conjecture  5.1  holds  t rue  for  open  networks  in 


f I ii  i I ihri  uni. 
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6 . Nonfeed b ack  N o n e x i t A r cs  in  Equilibrium  Open  J ack son  Networks 

In  this  section  we  deal  with  open  Jackson  networks  in  equilibrium; 
recall  (see  Sec.  1)  that  the  traffic  solution  6 satisfies  6 . =E|D . ( t , t+1 ] 1 

l l 

t>0,  for  any  i'-M.  Throughout  this  section  (a,b),  a^b,  designates  a 
traffic-nontrivial  nonexit  arc.  A major  simplification  will  be  attained 
in  the  impending  computations,  by  taking  advantage  of  the  identity 


(6.1)  <t>  (z  ,...,z  ,y)-6  (z  , . . . ,0. , . . . ,z  ,y)=$  (z  ,...,z  ,y)p.z. 

ol  m ol  j m ol  m 1 j 

which  follows  from  (2.3)  and  (2.8).  Furthermore,  the  equilibrium 

situation  brings  about  a substantial  reduction  of  combinatorics  bv 

enabling  us  to  collect  and  "complete"  terms  into  vanishing  expressions. 

Thus,  (5.2)  becomes  in  equilibrium 


(6.2)  D[<f>  (z  ,...,z  ,y)]=<£  (z  , . . . ,z  ,y)6  p ,z  (y-1). 

o 1 m o 1 m a ah  b 

A more  complex  application  of  the  "completion"  method  will  now  be 

n k 

exemplified  for  the  computation  of  Dip  (z, z ,v)(  z.  )|  when 

o 1 m k=]  lk 


u 11  ' k 

1 n-  m,  >1 , but  alJ  = { j , : h k^n ■ . Denoting  1 (z,,...,z  ,v) ( z . ) = 

- - k-  n ‘ k - - ol  m"k=llk 

and  with  the  aid  of  (6.2)  we  compute  (for  the  sake  of  clarity  we  shall 
explicitly  list  all  terms,  including  the  vanishing  ones) 


D [ f]  = 


}.  Y". . ( z . - 1 ) 
. i i 
i M 


+ !.  ff’.Z.i  q 


jrM-.I  1 1 1 1 j i J 

n n 

7.  . 


. (’  -D+)!  (’r-o)J.  q.('  -l) 


I 1 7.. 


+ ) \ Y(».z.o.p..(  ' -D+y  ) (f-D)  ,p..(  1 -1) 

i M jt  M-.l  1 ' 1 ''  7-j  i , M Jf'.l  1 " ' j 

i J*  i " 


f yP  z o p 

a a a nb  7. 


(v-1 ) 
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Next,  we  complete  the  right  side  of  the  above  to 
3 n &k 

K~r  <)>  (z  , ...,z  ,y)(  :i  z.  )=0  by  adding  and  subtracting  terms  of  the 
dt  ° 1 m ' k=i  Jk 

1 Zi 

form  YP.z.cf.q.( 1)  and  Yp.z.a.q.( 1)  for  ieM  and  jej  . We  obtain 

J J J J zj  J J J J z n 

after  some  algebraic  manipulation 
D[Y]  = 

ST  > • • • »zm*y)  ( S Z.  k) 
r)t  ° 1 m k=l  Jk 

+ I Y'*.q.(~ — l) (1-p.z.) 
jcj  J J zj  J J 

J n J 

z . 

+ I I Y°.P  . . (—  -1)  (1-P.z  . ) 

ieM  jcJ  ^ ^ zj 
n J 

+ ^aPab'b^"^  ' 

Finally,  after  collecting  terms,  the  end  product  of  the  calculation 
reduces  to 

n £, 

(6.3)  Dfi{i  (z  , ,z  , y)  ( i z.k)  ] = 

° 1 m k=l  Jk 

n p . m 

i 4>  (z  , . . . , z , y)  ( z . k) (1-p  . z . ) (a  .q  . (1-z. )+  ) i .p  (z  -z 

n m ii  7 i 1 i i i 1 ■ ill  i 


V J 


k=i  Jk  zj 


3 J J 


i=l  J Ii  i 1 

ifi 


)) 


+ <j>  (z  , ...,z  , y)  ( z.k)  <5  p z (y-1)  . 
o 1 m R=i  Jk.  a ab  b 


Let  /^(b)  be  the  set  of  all  paths  (b,j^,...,j  ) that  originate  from 

node  b and  have  length  n^l.  Let  / (b)=(it_M:  i(b,i, i ,,i)’  (b)  for 

- n 1 k - 1 k 

.some  l_k'n),  n 1 , be  the  n-neighborhood  of  b,  where  tor  n 0 we  define 
1/  (b)=lbl. 

O 

Finally,  let  ; ’ (b)  be  the  set  of  all  sums  whose  terms  have  the  form 
n 


(6.4)  i(/  z ,y)-C(y- 1 ) (z  z ,y)(  z ')  • (l-i  ./.)(/ 

I m o 1 m . , k , . k k e d 

k’  I k-  I 


where 


J 


1.  C is  a scalar 


2.  r>0,  r >0  are  integers 

3.  I.,I  (b) 

i 2 n 

4.  d,e£k  (b) 

n 

Lemma  6 . 1 


For  any  t(z  z ,y)f-,7  (b),  l^rr  L,  of  the  form  (6.4)  we  have 

i m n ~ 


D[t(z  ...  .,z  ,y)  ] . (b) 
1 m n+.l 


Proof 


Expand  ( -If  z,  ) FI  (1-p,  z ) into  a (finite)  power  series  in  the 
kelL  k kel2  k k 

z.  . Thus,  t(z, ,...,z  ,y)  is  a sum  of  terms  of  the  form 
K.  i m 

l. 

(6.5)  x'(z  ,...,z  ,y)=C'  (y-l)r(J)  (z  ,...,z  ,y)(  II  z )(z  -z  ) 

1 m oi  m kf  1 k e d 

for  some  C',  £ >0,  and  FI ,Ul  cy  (b).  Recalling  that  .liV  (b)  , we 
k-  1 2 n n 

shall  now  compute  D[ i ' (z, , . . . , z ,y)]  with  the  aid  of  (6.3);  for  the  sake 

1 m 

r A 

of  brevity  we  denote  C'(y-l)  <$>  (z  ,...,z  ,y)(  II  z k)  = y' 

oJ.  m . T k. 

kcl 

D[x' (z  , . . . ,z  , y) ]= 
i m 

1 m 
l y'z  — (1-P.z .) (O .q . (1-z. )+  ) O . p (z  — z . ) ) 
jeiljfe}  e zj  J J J 1 J i=l  J l1  1 i 

i^j 


+ Y’ze6aPabzb(y"1) 


II J fd  I 


y'z  — (1-f)  z . ) (o  q . (1-z  )+  j [a  .p  . (z-z.)) 

fl  zj  ) j i ) i i=l  .)  Ji  i j 


- Y,//aPabZb(y-|)  = 

1 H? 

}.  y'  - (l-P.z . ) (n  .q  . (1-z  )+  ),  o .p . . (z.-z, ) ) (z  -■/  ) 
jEl  Z.j  1 > ' 1 1 i = l 1 ” 1 J c (l 


\ Y1  (l-f'.z  . ) (cj  .q  . (1-z  . )+  ),  O p (z  -z  )) 
jcfe.dl  ’ 1 J 1 1 i-1  J 11  1 1 


+ Y’6aPabZb(y“1)(VZd)- 


An  inspection  of  the  right-most  side  above  will  verify  that 

D[t'(z  , — ,z  ,y)]eS  ,,(b).  The  Lemma  follows  from  the  linearity  of  D. 
i m n+i 


We  are  now  in  a position  to  state 
Theorem  6.1 

I.et  JN=(M,ra,a  ,P)  be  an  open  Jackson  network  in  equilibrium.  Then 
for  LriL 

•6.6)  Dn  [ t^o  (zi » • • • > zm»y)  ] = 

4>*(z  , . . .,z  ,y).z" 

. n o 1 m d 

a t 

,)  n_1 

+ 57  'I'Mz  , . . . ,y)  ) ( n (1-p  z )CJ  p -z.  ; 

')c  ° 1 m Pn_i(b)  k=!  Jk  >k  Jk  Vk+1  Jn  ln_i 


where  z ,y)=4>  (z  , . . . ,z  )e 

t 1 m o 1 m 


^Pabt(y_1) 


, and  R „c5  „ (b) ; the 

n-z  n-Z 


sums  V here  and  elsewhere  are  taken  over  all  (j  , . . . , i )eP  .(b), 
P ,(b)  lnn-1 

n-i 

unless  otherwise  specified. 


From  (6.2) 

(6.7)  D[<Mz1,...,zn],y)]=  ^ 4>*(Zj Vy)'Zb  * 

The  proof  will  proceed  hy  induction  on  n=2,...,L. 

For  n=2  we  compute  from  (6.7)  with  the  aid  of  (6.3) 

D2[VZl’-**’Vy)1=D[37  ^zl Vy)zbl  = 

iTt  K(Z1 Zm*y)(1-PbZ!,)(,,bf’b(1-Zb)+.^',|Ai(zi-7|,)) 

i =1 


+ Jt  ^(zr--*’Zm’y)zhVnbZb(y-,)  = 


i 
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r ...  . 2 

— = 4>*(z  , . . . ,z  ,y)z 
l o 1 m n 

3 t 


+ •^*5^! z~’y>  ( « n-(»,  z.  >r.  p. 


/ x (b)  k=l 


k 'k  'k  Vk+1 


)(z.  -z.  ) 


l2  ) j 


+ R 


for  some  R zS  (b), which  establishes  the  induction  basis, 
o o 


Now,  assume  that  Dn[<t  (z,  ,...,z  ,y)  ] has  the  requisite  form  (6.6) 

o 1 m 


for  2^  n^  L , and  apply  D to  D [4>o  (z^  , . . . , z^,y)  ] . From  (6.3)  the  first 


expression  in  (6.6)  Is 
n 

(6.8)  I)[  <|>*(z z ,y)z"]  = 

n o 1 m b 


<)  t 

v n+1 


7tnTT  K(zi Vy)7b+1 


+ R 


(1) 


for  some  (b)‘  !J  . (b) 

l 1 n- 1 


Consider  the  second  expression  in  (6.6).  By  the  induction  hypo- 
hesis,  each  term  T . . . . in  the  sum  associated  with  the  path 

(J1’---’V 


(j  )eP  , (b)  belongs  to  S (b) . For  each  fixed  (j  )rp  .(b), 

I n n-l  n I n n-1 

n-1 


expand  FI  (1-p.  z.  ) in  the  corresponding  T(  , , as  in  Lemma  6.1, 

k=l  'k  ’k  ''•11 V 


and  set  e= j and  d=j  Thus  for  each  term  T 1 of  the  form  (6.5) 

' n n-i 


such  that  i,.  . . =)l'  (we  use  the  results  and  notation  ol  l.cmnvi  6.1) 


l> ( T * (Zj  , . . . ,zm,y)  J = 


m 


Y'd-C  z )l  a p j(z.  -z  ) 
JnJni=i  JnJn  A Jn 


ifi. 


+ R' 
n- 1 


lor  some  R'  >.!!  .(b).  Note  that 

n-1  n-1 


n-1 


/.  Y ' = 


|)*(z.,...,z  ,y)  n (1-p  z.  )'!.p.  . the  sum  being  taken 


-It  'o'  1 


k=1 


'k  'k  1 'k'k+l 


over  all  i'  whose  sum  comprises  I.  Tims,  smnmiri)',  over  all  such  i * , 


, I )'  /'  . d>) 

n n- 1 


J i 


we  oh  lain  lor  each  ( j 


* • • 
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D[T(3l J„)<zi v^1 


if>*(z  ,...,z  ,v)(  n (1-p  z )0  p . Kl-P,  z ) I P,  p.  i^zi-71  1 
3t  o 1 m R=1  Jk  Jk  JR  Jk  lk^L  >n  Jn  i=l  3n  Jn  n 

+ R 


n-1 


m n 

~~  <p*(z  ...  .,z  ,y)  l ( li  (1-P..  7 )fT  p )(z  -7  ) 

r>t  ° 1 i ,=1  k=]  ’k  ]k  'k  'k  1 k+1  'n+1  Jn 

n+1 


4-  R"  . 
n-1 


for  some  R . i u . (b) . 

n-1  n-1 

Summing  again  the  above  over  all  wt‘ 


n-1 


(6.9)  Dk- 


- <j)*(z  , . . . ,z  ,y)  l ( II  (1-P  z 

t ° 1 m P . (b)  k=l  k 1 

n-i 


)0  . p . . ) (z . -z.  ) 1 

k 'k  'k’k+l  'll  ’n-1 


4z  0*(z ,,...,z  ,y)  l ( n (1-P  z )o  p )(z  -z  ) 

at  0 1 m P (b)  k=l  Jk  ]k  Jk  ]kJk+l  n+1  ’n 

(2)  n 


+ R . 
n-1 

(2) 

fo  r so  mi'  l<  . 1 .0>). 

n-1  n-1 


Finally,  from  Lemma  6. 


(6.10)  1)  | R ' ,(b). 

n-2  n-i 


The  induction  step  follows  from  (6.8),  (6.9)  and  (6.10)  by  linearity  of  I) 


The  main  result  now  follows. 

Theorem  6.2 

l.et  .IN-  (M, o,o  , P)  he  an  open  Jackson  network  in  ecpi i I i hr  i um.  Then 

iK  . (t)l  is  not  a Poisson  process, 
an  to 

Proof 


For  n=L,  set  in  (6.6)  z^  = l for  all  i?<a.  Since  aeVj  (b)-Fj  -,(6), 
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Rj  ^ vanishes  along  with  all  terras  corresponding  to  paths  (j  ,...,j  ) 1 (6) 

with  j ^a,  leaving  only  terms  that  correspond  to  all  minimal  cycles  of 

L 

(a,b) . Thus, 

(6-U)  ^VVy)  = 

<},o(za’y)('aPab(y-1))L 


+ ).  ( IT  (l-p  )a  , p ) (z  -1)  , 

C . (a,b)  k=l  'k  'k  ]k'k+l  a 
min 

wll(‘rt'  ).  is  taken  over  all  (a,  j, j , ) r C (a,b). 

C . (a,  b)  1 '•  min 


On  setting  z =0,y=0  above, we  obtain 

3 

L 

(6.12)  ~~  P[  B (0)=0,K  , (0)=0]  = 

0 tL  d ab 

P[Ba(0)-0.Kab(°)-01(-«aPal))1 

L-l 

l ( II  (l-p.  )o.  p.  . ) 

C . (a,b)  k=l  'k  *k  k k+1 
mm 

Now,  if  fK^^Ct) } q is  a Poisson  process,  then  from  Corollary  3.2 


(6.13) 


, PfBa(0)=0,Kab(0)=0]= 
lim  P f 15  (0)=0]e''aP;ll),(-.';p  )'  = 


,,[Ba(0)=0’Kah(0)=nl(_<Sapab) 


But  for  each  (a,  j j , . . . , jj  )(  i’min(n,b)  we  have 


II  (1-0 . )'J,  P.  . 0. 

k=  1 'k  'k  'k'k+l 


As  (’  . (n,b)/(|>  by  our  assumption,  (6.12)  contradicts  (6.13)  whence  the 
m I n 


Thi’O  rein  loll  ows  . 


7. 


Conclusions 


The  intuitive  essence  of  the  preceding  results  may  be  loosely 

described  as  follows.  The  traffic  process  on  arc  (a,b)  is  statistically 

dependent  on  all  state  indicators  of  those  nodes  which  it  affects  (all 

nodes  accessible  from  node  b) . These  dependencies  radiate  about  node  b 

in  a sort  of  a ripple  effect;  they  manifest  themselves  for  nodes  which 

.^n+1 

are  n-distant  from  b through  —r  <p  (z z ,y).  The  existence  of  a 

..  n+1  o 1 m 

d t 

minimal  cycle  for  (a,b)  enables  these  dependencies  to  propagate  to  node 
a,  thereby  leading  to  a contradiction  to  Theorem  3.1. 

As  a consequence  of  Theorem  6.2,  one  may  summarize  the  foregoing 
discussion  concerning  traffic  processes  in  equilibrium  open  Jackson 
network,  as  follows. 

Conclusion  7.1 

Let  JN=(M,a,CT,P)  be  an  open  Jackson  network  in  equil ibrium, and  let 
(a,b)  be  a traffic-nontrivial  arc.  Then  the  following  are  equivalent 
statements: 

a)  ( a , b ) is  an  exit  arc. 

b)  IK  . (t)l  ,,  is  a Poisson  process. 

nb  t-'O 

c)  IK  . ( t ))  is  a renewal  process. 

ah  1 O 

<l)  For  each  t 0,  K , (t)  is  Poisson  distributed  with  parameter  p . t. 

ab  a ah 

e)  For  eacli  t>0,  B (t)  and  K . (t)  are  independent  random  variables. 

— cl  <1  O 

I)  There  is  N' M,  a N,  such  that  for  eacli  t'O,  K . (t)  is  independent 

ab 

of  I ( ) . ( t ) : i r N ) ; in  fact,  N - i : i . - > a ; . 

Another  conclusion  reflecting,  on  Theorem  tell:,  us  about  the 

limitations  inherent  in  equilibrium  dccomposit  ions  ol  open  Jackson  nit 
works . 
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Conclusion  7.2 

The  decomposition  in  Theorem  2.2  of  an  open  Jackson  network  in 
equilibrium  into  Jackson  network  components  is  maximal  in  the  sense  that 
no  refinement  of  the  underlying  partition  wil]  yield  Jackson  subnetworks 
( i n equ  i I i hr  i um)  . 

1 1 

Next,  let  us  extend  the  notion  of  an  exit  arc  to  any  superposition 
of  traffic  processes  on  arcs  emanating  from  the  same  node,  as  follows: 
the  superposition  above  is  said  to  have  the  exit  property, if  all  its 
constituent  arcs  are  exit  arcs;  otherwise,  it  has  the  nonexit  proper-  ty  • 

In  the  latter  case,  the  minimal  cycle  associated  with  it  is  the  smallest 
over  all  nonexit  arcs  in  the  superposition.  We  point  out  that  all  the 
foregoing  results  will  still  hold  nrutatir,  rrrutnndir,  for  such  superpositions. 

Going  back  to  Conjecture  5.1  and  comparing  it  with  Corollary  3.2, 
we  see  that  the  sum  d = !)  d(a, j , . . . , j.  ) | has  a useful 

C . ( a , b ) 

min 

interpretation;  it  may  be  construed  as  a heuristic  deviation  of 

fk  . (t)}  from  the  Poisson.  Moreover,  each  term  in  the  sum  mav  be 
ab  t 0 

interpreted  as  the  deviation  due  to  the  respective  minimal  cycle  of  (a,b). 

In  particular,  for  open  Jackson  networks  in  equilibrium,  this 
deviation  takes  on  the  form 


To  support  this  view  we  point  out  that  in  view  ol  l> ) in  Conclusion 

/ . I , K ( I ) ,,  is  a Poisson  process  ill  (’  . (a , l> ) , and  con  sc- 
ab I _(•  min 

qircnt  I v ill  d (>. 


Equation  (7.1)  provides  some  clues  to  the  behavior  of  the  deviation 
as  a function  of  the  network's  structural  parameters.  A number  of  quick 
observations  now  follow: 


1.  The  deviation  magnitude  is  affected  by  network  connectivity  in  the 

sense  that  it  increases  as  the  number  of  minimal  cycles  increases. 

7.  The  devial ion  component  due  lo  a minimal  cycle  is  aflecled  by 

distance, in  the  sense  that  it  decreases  as  the  minimal  cycle  length 

I.  increases,  provided  the  service  efforts  . p along  ares  (i  i ) 

'k'k'k+l  k kH 

are  less  than  one. 

3.  The  deviation  component  due  to  a minimal  cycle  is  affected  by  tin- 

cycle's  "strength".  It  decreases  as  the  cycle  gets  "weaker", in 

the  sense  that  the  switching  probabilities  p.  . decrease  in 

]kJk+l 

magnitude . 

4.  The  deviation  component  due  to  minimal  cycle  are  affected  by  the 

work  load.  It  decreases  as  the  nodes  tracing  it  £et  busier,  i.e. 

as  the  traffic  intensities  o approach  one. 

'k 

As  l matter  ol  fact  on  approaching  saturation,  tin  respei  live 

deviations  tend  to  0.  Consequently,  under  heavy  traffic  conditions, 

over  all  paths  in  C . (a,b),  !K  . (t):  _ is  approximately  a I’oisson 

nan  ah  t 0 

p roi  i ss . 

Finally,  we  ment  ion  that  analogous  versions  ol  the  results 
contained  in  this  paper  can  he  conceivably  derived  lot  a general i zed 

notion  ol  trallie  processes  in  Markovian  systems.  The  tcasihilily 
ol  smh  generalizations  is  currently  under  investigation. 


II  the  point  " is  added  to  t he  state  space  ol  each  node,  then  the 
equilibrium  notion  will  include  the  case  P|Q,-,  (t  )=- ■ ) 1 (saturation). 
Not  ice  that  a saturated  node  is  indist inguishable  Irom  the  environ- 


ment  source  because  all  customer  streams  departing  that  node  are 
mutually  independent  Poisson  processes  which  are  also  independent  of 
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